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Abstract Halo coronal mass ejections (CMEs) were found to be significantly faster than 
normal CMEs, which was a long-standing puzzle. In order to solve the puzzle, we first 
investigate the observed properties of 31 limb CMEs that display clearly loop-shaped 
frontal loops. The observational results show a strong tendency that slower CMEs are 
weaker in the white-light intensity. Then, we perform a Monte Carlo simulation of 20 000 
artificial limb CMEs that have average velocity of ~523 km s^^. The Thomson scattering 
of these events is calculated when they are assumed to be observed as limb and halo 
events, respectively. It is found that the white-light intensity of many slow CMEs becomes 
remarkably reduced as they turn from being viewed as a limb event to as a halo event. 
When the intensity is below the background solar wind fluctuation, it is assumed that they 
would be missed by coronagraphs. The average velocity of "detectable" halo CMEs is 
■~922 km s~^, very close to the observed value. It also indicates that wider events are 
more likely to be recorded. The results soundly suggest that the higher average velocity 
of halo CMEs is due to that a majority of slow events and a part of narrow fast events 
carrying less material are so faint that they are blended with the solar wind fluctuations, 
and therefore are not observed. 

Key words: Sun: coronal mass ejections (CMEs) — Sun: activity — methods: numerical 
— solar-terrestrial relations 



1 INTRODUCTION 

Since th e first r e cogni tion of coronal mass ejections (CMEs, initially called coronal transients) by 
OSO-7 dTousevl Il973h . more than 10000 such energetic events have been identified by ground- 
based and space-borne coronagraphs. Althoug h remarkable progress h ad been made before, the Large 
Angle Spectrometric Coronagraph (LASCO) jBrueckner et al.U 19951) . which is aboard the Solar and 



Heliospheric Observatory (SOHO) spacecraft ( Domingo et al.l 1995b launched at the end of 1995 



revolutionized our understanding of this eruptive activity for its large field of view (FOV), increased 
sensitivity and dynamic range. CMEs are often associated with solar flares and filament eruptions 



Chen et all 120061; IChen & Zond. l2009h . leading to large-scale coronal disturbance s like EIT waves 



ChenL l2006i 120091) 7 and even triggering a sympathetic CME (ICheng et all 12005 ah . A typical CME 



exhibits the 3-part morphology: a frontal loop, which is foUowed by a dark cavity with an embed- 
ded bright core (e.g..lllling & Hundhauseniri985HDere et al.Lll997l) . It expels approximately lO^'^-lO^*' 



gram (e.g.. lWebb et al.l 19951) of plasma into interplanetary spa ce with velocity ranging from tens to 



300 km s-^ (e.g..lSt Cvr et all 1 19921; iHundhausen et al.L 1 1994 . and with kinetic energy up to ^10 
erg jBurkepile et alT 2004 ). At the same time, magnetic field and its intrinsic magnetic helicity are 



ejected into interplanetary space, which plays an essential part in completing the global magnetic field 
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reversal between successive solar cycles (e.g.. lZhang & Lowl 2005h. Thev can po tentially give rise to 
hazardous terrestrial effec ts, such as solar energetic particles (Reames et all [19991) . type II radio burst 
jGopalswamv et al"ll2009l) . geomagnetic storms (e.g. . iGosling et alT 1990l) . ionosphere disturbance, and 
polar aurorae. The information of the magnetic field, mass, and velocity o f CMEs is very important since 
they determine the geomagnetic effectiveness dGopalswamv et allbOOTb . 

As a special type, those surrounding the occu lting disk, i.e., with an apparent angular width of 360°, 
are called full halo CMEs dHoward et all Il982l) . Compared with normal events with apparent widths 
between 20° and 120°, they are generally believed to be nothing special except that they propagate in a 
direction close to the Sun-Earth line, either toward or away from the Earth. However, it has been noticed 
that th e average app arent velocity of halo CMEs is fairly higher than that of normal CMEs (Webb et all 
I1999I) . For example. lYashiro et al.l (|2004|) compared ~7000 events in the period from 1996 to 2002, and 
found that the ave rage app arent velocity of halo CMEs is twice larger than that of normal CMEs. Such 
a difference made lLara et al. (2006) propose that halo CMEs are of a special type. Realizing that white- 
light emission of CMEs comes mainly from the Thomson-scattering of photospheric radiation, which is 
much weaker for plas ma near the Sun -Earth line than that near the plane of the sky at the same projected 
heliocentric distance, lAndrewsl(l2002l) suggested that there exists a mass cut-off, above which CMEs are 
bright enough to be detected, while many dim, presumably slow events, are missed by coronagraphs. 
This would make the average apparent velocity of halo CMEs much higher than that of the normal type. 

Suc h a conjecture can be validated if slower CMEs are systematically fainter in brightness. For this 
purpose. lcheng et"ai1 ( l2005bl) studied the relationship between the apparent velocity and the white-light 
intensity for halo CMEs. The two parameters did show a marked positive correlation, which provides 
an indirect support for the mass cut-off conjecture. Compared to limb events, halo CMEs should travel 
a longer distance to be observed in the FOV of coronagraphs, which has two effects making halo CMEs 
significantly fainter in white light. First, the intensity of the incident emission from the photosphere is 
lower Second, the number density of the CME front be comes smaller. Despite that the cross-section 
of the Thomson scattering of halo CMEs increases (see iBiUingslI 19661) . their white-light emission is 
reduced greatly compared to the limb CMEs events that are observed at the same projected distance in 
the plane-of-the-sky. Therefore, it is expected that many halo CMEs, especially the slower events, could 
be so faint that are missed by coronagraphs. In this paper, we collect a sample of 31 limb CMEs free 
from projection effect during the SOHO Mission to confirm the velocity-brightness relation. A Monte 
Carlo simulation is further performed to quantitatively testify whether the observed high value of the 
average velocity of halo CMEs is due to that many slow events are missed by coronagraphs. 

This paper is organized as follows. Data analysis and the results are presented in Section |2] The 
Monte Carlo simulation and its result are shown in Section[3] followed by the discussion on the projec- 
tion effects in Section|4] We summarize the results in Section|5] along with some discussions. 



2 DATA ANALYSIS AND RESULTS 



Located at the inner Lagrangian point (Li), SOHO has been monitoring the vigorous Sun for 14 years. 
Three coronagraphs, CI, C2, and C3, which constitute the LASCO instrument, have FOVs of l.l-3i?0, 
2-6i?0, and 4-32i?Q, respectively, where Rq is the solar radius. The routine observations provide a 
huge database for CME research. 

In or der to select limb events, we carefully checked the movies that are composed of LASCO/C2 
and FIT jDelaboudiniere et all Il995h difference images from the NASA CME catalog between 1997 
January and 2005 December The events with associated flares or filament eruptions occurring beyond 
the longitude of 50° are collected. Moreover, only the CMEs that have definitely loop-shaped leading 
edge are considered. It should be kept in mind that the bright loop is not a simple tube but the projection 
of the dense front of three-dimensional bubble-like structure. To select events with clear loop tops, those 
CME events with their legs being much brighter than the tops are ignored. Besides, some events whose 
white-light intensity increases with height or can poorly be fitted with a power-law function with height 



http://cdaw.gsfc.nasa.gov/CME_list 
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Table 1 List of basic properties of the 31 loop-shaped limb CMEs, including 
the date, time of first appearance in C2 FOV, central position angle (CPA), 
angular width (AW), and the linear velocity (V). 



Date 


Time 


CPA 


AW 


V 




(UT) 


(Deg) 


(Deg) 


(kms^^) 




1(5. jU 


/o 


OQ 

yy 


Dvyy 


1 QQ7/1 /^'^ 


U(5.Z / 


Z /J 


QQ 
yo 


414 




14. JO 


Zoo 


/4 


Z4u 


Lyyo/\jZ/Zj 




7/1 

/4 


OJ 


Zoy 


1 QQfi/1 1 /no 
Lyyo/ i l/\Jy 


Ol -^A 


1 

10 


74 


AAA 


1 QQO/n'^/1 7 


nn.cn 
UU.jU 


Zyj 


1 Ij 


^o^ 
jyjj 


zUUU/U(5/ iZ 


1 J.j4 


ZD4 


11/ 


AQQ 

4yy 


ZUUU/Uo/ZZ 


Zj.UD 


1 7Q 

1 /y 


jy 


All 

4J1 


OOOO/I 1 /II 
ZUUU/ 1 1/Z / 


Zj.j4 


1 7T 
IZj 


<7 
J / 


4/4 


0OO1 /Ol /TA 
ZUUi/Ul/Zu 


10. UO 


cc 

J J 


111 
111 


oycS 


0OO1 /O'^/TQ 
ZUUi/Uj/Zo 


O^-^O 
Zo.jU 


yo 


A 1 
41 


<A0 
j4Z 


ZUUi/UD/ i J 


OO-OA 


77Q 

Z /y 


oz 


44/ 


0OO1 /OC/QO 

ZUUl/UcS/JU 


nQ.cn 


1 7Q 

izy 


(5U 


40Z 


0OO1 /0Q/0< 

zuui/uy/uj 


10. uo 


ZjZ 


1 07 


<1C 
JJcS 


2002/01/10 


00:30 


236 


61 


377 


2002/03/12 


23:54 


112 


82 


535 


2002/03/15 


02:06 


64 


65 


272 


2002/04/18 


06:26 


162 


64 


552 


2002/09/18 


14:54 


279 


99 


512 


2002/11/04 


12:30 


17 


114 


509 


2002/11/08 


11:30 


298 


69 


424 


2003/01/03 


11:30 


283 


88 


521 


2003/01/20 


18:30 


315 


105 


733 


2003/03/21 


10:54 


54 


66 


481 


2003/12/08 


13:31 


228 


68 


464 


2004/05/03 


00:50 


113 


112 


464 


2004/07/10 


13:54 


270 


78 


477 


2004/08/18 


17:54 


258 


120 


602 


2004/08/27 


09:30 


261 


70 


554 


2005/03/14 


08:00 


259 


105 


849 


2005/09/04 


14:48 


286 


86 


1 179 



are also excluded (possibly due to that they are undergoing acceleration). As a result, the sample finally 
consists of 3 1 well-defined events, whose image quality indices are > 4 (labeled as "good" events in 
the CME catalog). To our understanding, the events with high quality appear less diffuse and have 
sharp contrast to the background corona so that the height-time measurement is more precise. The basic 
properties, including the observation date, time of first appearance in the C2 FOV, central position angle 
(CPA), angular width (AW), and the linear velocity (V), are listed in Table [T] 

Pre-processing of the LASCO/C2 data, such as correction of flat field and removal of dark current, 
is conducted by using the standard program cljoalibmte.pro in the Solar SoftWare (SSW). For each 
event, we take the white-light image just before the first appearance in the C2 FOV as base image to get 
base difference intensity of the CME frontal loop in the ensuing three snapshots. 

It is noted that the observed white-light intensity is the emission integrated along the line-of-sight 
(LOS), which represents the total mass of unit area in the plane-of-the-sky, i.e., column density. We call it 
integral intensity. In order to better characterize volume density enhancement of the CME leading loop, 
we define another parameter - normalized intensity, which is the integral intensity divided by the LOS 
depth in unit of Rq . Since the 3D topology of CMEs is unclear from a single viewing direction, several 
versions of the cone model have been proposed based on the fact that the app arent angular widths of 
many CMEs reniain almost constant during their propagation i n the corona (e.g^ Anzer & Polandl 19791 : 
iFisher & Mu"nrolll984l;lMichalek et all 120031: IXie et alll2004 . According to lSchwenn et all (l2005h . flie 
projected geometry that can well reproduce the kinematic properties of CMEs with the cone angle 
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Fig. 1 A sketch of the plane-of-the-sky-projected axisymmetric cone model applied to the 
analysis of loop-shaped limb CMEs at their early stage of propagation in the corona. The 
cone OMTNO outlines the CME leading edge, while 7 and ji are the central position angle 
and the angular half-width, respectively. Line CY perpendicular to OT equals to half LOS 
depth of Y. 

between 40° and 80° is displayed in Fig. [T] The frontal loop is concentric with the solar disk. The 
position angle of cone axis OT is denoted by 7, and the angular half-width by /3, respectively. The LOS 
depth of white-light emission, which is an unknown parameter, is supposed to be of the same size as the 
transverse expansion according to the cone model. In the right triangle /\OYC {Y denotes a point along 
the hne 07 within the frontal loop), OC = OT and YC = ^/{OCY - {OY f = ^{OT f - {OY f. 
Under the axisymmetry assumption of the CME 3D morphology, the LOS depth at Y equals to 2YC. 

For each CME, the normalized intensity of the loop-top at three moments is obtained and fitted 
with a power-law function /„ ~ r~^, where r and k are heliocentric distance and the power index, 
respectively. As an example, we plot the variation of the normalized intensity /„ with r for the 2000 
August 22 CME event in Fig. ID which corresponds to fc = 3.60. In the same way, we derive k for 
each loop-like event based on near-perfect power-law regression. The distribution of the index, which 
is averaged at 3.85, is shown in Fig. |3] i.e., on avarage, /„ decreases with the heliocentric distance 
r as In ~ 7--3-85. As mentioned in Section [T] two factors contribute to the decrease of /„ with r. 
One is the incident emission from the photosphere, the other is the number density of the CME front. 
It is known that the incident emission decreases with r as r~'^, therefore, it is derived from the 31 
limb events that the number density of the CME frontal loop decreases with height as ^ ^.-i 85^ ^ 
little more slowly than the first factor Since a flare- or filament-related CME usually undergoes three 
phases: initiatio n phase, fast acceleration phase in the inner corona (< 3Rq), and propagation phase with 
constant speed (IZhang et al.ll200lll2004l) . the white-light intensity interpolated at 3Rq, both before and 
after normalization, is calculated to check the relation between CME velocity and brightness. For the 
2000 August 22 event, the integral and normalized intensity at 3Rq is 1.70 x 10^^ msb and 4.47 x 10~^° 
msb, respectively. Here, "msb" is short for "mean solar brightness". 

The relation between CME velocity and integral intensity at the frontal loop center (which is also 
the brightest point) is shown in the left panel of Fig. ID A clear tendency is seen that the integral intensity 
of the CME leading loop increases with the CME velocity. The correlation coefficient is as high as 0.62. 
Apparently, the 3 1 limb CMEs all fall in the domain bracketed by the two parallel soUd Unes in the panel. 
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Fig. 2 Variation of the normalized intensity /„ of the CME leading edge with the heliocentric 
distance r for the 2000 August 22 event. The variation is fitted with a power-law profile 



where k = 3.60. 
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Fig. 3 Scatter plot of the power index k in the 31 loop-like events. The average value of k is 
3.85. 



The relation between CME velocity and the normalized intensity at the frontal loop center is shown in 
the right panel of Fig.|4] The correlation coefficient increases to 0.72 after normalization. Moreover, the 
scattering of the data points is reduced compared with the left panel. 

The positive c orrelation between th e white-light intensity and the velocity of the limb events con- 
firms the result of ICheng et al.l (l2005bl) . and provides indirect evidence in favor of the conjecture of 
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Fig. 4 Scatter diagrams showing the relationship between the velocity and the integral {Left 
panel) and normalized {Right panel) intensity at the loop center of the 31 sample events. The 
solid lines bordering the data points in the left panel are parallel to the fitted line from these 
points. The correlation coefficients are displayed at the upper left of the panels. 



I Andrewsl (l2002b . who proposed that the high average halo CME velocity is due to that some slow events 
are neglected by coronagraphs. To quantitatively justify this viewpoint, we carry out a Monte Carlo 
simulation in the next session on the basis of the above correlation. 

3 MONTE CARLO SIMULATION 

As mentioned in Section[T] the white-light emission scattered by the CME front in the halo events would 
be greatly reduced, compared to limb events. Once the white-light intensity falls below the fluctuation 
level of the background solar wind, it is likely that a halo CME event is missed by coronagraphs. In this 
section, using Monte Carlo simulation of 20 000 artificial CMEs, we attempt to quantitatively estimate 
how these effects change the average velocity of the observed halo CMEs. 

Figure |5] depicts the geometry of a limb {left panel) and a halo {right panel) CME, based on which 
we now check how a limb event changes in the apparent velocity and the white-light intensity when it 
is observed as a halo event. From the geometry, it is easy to see that the apparent velocity decreases by 
a factor of sin (3 when a limb CME is observed face-on as a halo event, where (3 is the true angular half- 
width. Line PQ in each panel stands for the longest LOS depth of the cone shell. The integral intensity 
in the limb and halo cases are labeled with II and Ih, respectively. 

The change in brightness caused by Thomson-scattering effect is a little complicated. For compa- 
rability, the white-light intensity is calculated at a projected distance of ?>Rq for both limb and halo 
events. Note that the real heliocentric distance of halo CME become Si?©/ sin/?. The frontal loop is 
assumed to have a radial width of O . Si?© at 3i?0 and increases by 1/ sin /?. The integral formula for the 
LOS white-light intensity jBillingsl 1 19661) is then expressed as 

1 r'^ 

I ^ -naoJoR n[{l - u){2C cos^^ e - A) +ui2Dcos-^ e - B)]de, (1) 

2 Jp 

where the density n is equal to for the limb case and to ri3(l/ sin/3)~^-^^ for the halo case, is 
the plasma density of the CME front at r = 3Rq, A = cosfisin^ fl, B ^ -0.125[1 - 3sin^ - 
+ 3siii^»)ln ^+7^" ], C = (4 - 3cosl7 - cos^n)/3, D = 0.125[5 + sin^ r2 - ^^^(5 - 
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To the observer To the observer 

Fig. 5 Two sketches illustrating how white-light coronagraphs view CME from an edge-on 
{Left panel) and face-on {Right panel) perspective. 

sin'^ ri) hi ii^Y^], sin57 = Rq/t, u ~ 0.6 is the limb darkening coefficient, ctq is the Thomson- 
scattering cross section, Jq the photospheric radiation at the solar surface, and the angular distance 
between the plane of the sky and the line connecting the Sun center with any point along PQ. Besides, 
we assume that the plasma density is uniform in the frontal loop for simplicity. For different angular 
half-width p, the ratio of the integral intensities of the halo and limb CMEs, Ih /II, is plotted in Fig. 
|6l It is seen that for a limb CME, the white-light intensity is reduced by 1-2 orders of magnitude as it is 
observed as a halo event. From the figure, we can also infer that the wid er a halo CME rea lly extends, the 
higher likelihood it would be detected. Such a result is consistent with lFainshteinI ( |2006|) . who reported 
that the real angular widths of halo CMEs with an average value of >60° are relatively larger than that 
of normal CMEs. 

In the next step, we perform a Monte Carlo simulation to generate a sample of 20 000 artificial limb 
CMEs that follow the real velocity and angular width distributions. Generally, both the CME velocity 
and the angular width suffer from the projection effects. In order to eliminate the projection effects, 
iBurkepile et alj ( |2004|) studied the statistical properties of limb CMEs, whose observed velocity and 
angular width are almost real. Using their data set, we display th e velocity d i stribut ion of the 72 events 
with velocity > 100 km s^^ in the left panel of Fig.|2l Similar to lAoki et alj (l2003h . the profile is fitted 
with a log-normal function, which results in fv = 0-19 cxp[— (In V — 5.76)^/0.48] + 0.01 as displayed 
by the dashed line. The angular width distribution of the 105 events with /? > 10° is illustrated in 
the right panel of Fig. [T] The histogram is fitted with a Gaussian function = 0.14exp[— (2/3 — 
41.95)2/7.02] + 0.03, as shown by the dashed line. The 20000 artificial CMEs are assumed to follow 
these distributions, with the average velocity being 523 km s~^ and the average angular width being 
53.8°, respectively. 

Then, for each velocity interval, data points with different angular widths are randomly distributed 
within the two bounding lines shown in the left panel of Fig. ID The 20000 artificial limb CMEs are 
scatter-plotted as red crosses in Fig. [8] in the case of k ~ 3.85. Note that only the randomly selected 
10% of the data points are shown in order to make the diagram clear. The thick solid horizontal line 
marks the 3(7 level of the white-light noise of the background solar corona at 3Rq, below which a 
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Fig. 6 The ratio of integral white-light intensity of halo CME to that of normal CME as a 
function of /3 at a projected heliocentric distance 3Rq. 
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Fig. 7 Velocity (Left panel) and angular half-width (Ri^ht panel) dis t ributio ns of limb CMEs. 
The histograms are from the observational results of lBurkepile et al.l ( |2004|) . while the dashed 
lines are out fitted profiles, a log-normal function for the velocity and a Gaussian function for 
the angular half-width. 



data point is considered as unobservable. Here a — 2.33 x lO^^'^ msb is the standard deviation of 
the background fluctuations in the LASCO/C2 images at 3Rq. For each data point, as the event is 
observed as a halo CME, its apparent (plane-of-the-sky) velocity decreases by a factor of sin /?, while 
its white-light intensity drops by Ih/Il (see Fig. |6]l. After such corrections, the new data points of 
the corresponding halo CMEs are plotted as blue diamonds. It is found that if the artificial events are 
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Fig. 8 Scatter diagrams of the randomly selected one tenth of the 20 000 artificial events. Red 
crosses and blue diamonds represent limb and halo events. The horizontal Une marks the 3cr 
level of the background solar wind fluctuations at 3Rq. 



Table 2 Parameter survey of the average velocity of halo CMEs {V, in unit of km s 
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observed as halo CMEs, the white-light intensity of a majority of the sample falls below the 3cr level, 
and therefore the corresponding CMEs are considered to be "missed" by the LASCO coronagraph. The 
average velocity (V) of the "visible" halo CMEs is calculated to be 922 km s^^. 

There are several assumptions about the properties of the CME frontal loop in the simulation, e.g., 
the radial width (h), the density variation with height (r"*^'^^^-'), and so on, where the density variation 
with height was derived from the 31 well-defined CMEs. To show how the result is affected by the 
assumptions, we performed a series of simulations with different h and k. The corresponding V is 
displayed in Table|2] with the unit of km s~^. 

4 PROJECTION EFFECT 

As one of the puzzles in CME research, it is recognized that halo CMEs are much faster than normal 
events, although in principle, the difference between halo and limb events is only the direction of prop- 
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Fig. 9 The observing geometry of halo CMEs in the front-side {left panel) and back-side 
{right panel) cases. AB signifies the plane of the sky, and S the observing spacecraft (SOHO). 




Fig. 10 The relationship between V^bs/K-eai and /3 for fore-side {Left panel) and back-side 
{Right panel) halo CMEs. Different line styles (solid, dotted, dashed, and dash-dotted) repre- 
sent the relationship at different projected heliocentric distances (3, 6, 15, and 'HRq). 



agation. Here we discuss whether the simple projection effect is related to the high velocity puzzle of 
halo CMEs. 

Since halo CMEs propagate either toward or away from the Earth as illustrated in Fig.|9] the leading 
edge (e.g., point D) is projected to the plane of the sky (i.e., point E) to calculate the propagation ve- 
locity of the CME. The left/right panel of Fig. |9]depicts the propagation geometry for a front-side/back- 
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side CME with real angular half-width j3. The real heliocentric distance of the leading loop, r, is equal 
to the length of OD, while the apparent heliocentric distance, x, is equal to the length of OB. The 
SOHO spacecraft, which is located at Li with a heliocentric distance of L = 212i?0, is marked with 
S. According to the geometrical relation, the ratio between the apparent velocity (Vobs) and the real 
velocity (14eai) as a function of /3 and x is expressed as 

V IV _<^i<^ _ \ [l + x/{Lta.TLl3)fs\nP, front-side, 
V'obs/ Kroal - / - I [1 _ :r/(Ltan/3)]2 sin/3, back-side. 

The relationship between the velocity ratio and (3 in both cases is displayed in Fig.[TO](left panel for 
the front-side events and right panel for the back-side), where different lines correspond to the relation 
at different apparent heliocentric distances within LASCO FOV. It is seen that for both types of events, 
the apparent velocity is always smaller than the real one within the LASCO FOV. For the front-side 
halo CMEs, the velocity-reducing factor increases from 0.2 to 0.88 as (3 increases from 10° to 60° when 
the projected heUocentric distance is 3i?0. At a larger distance, the reducing factor increases. For the 
back-side events, the reducing factor is slightly smaller than that in the front-side case at 3i?0. However, 
the factor decreases with distance as indicated in the right panel. 

The above theoretical analysis reveals that the projection effect results in a smaller apparent velocity 
of halo CMEs. Therefore, the pr ojection effect cannot resol ve the high ve l ocity p uzzle of halo CMEs. 
Such a result is consistent with iMichalek et all ( l2003h and [Howard et all ( |2008|) . who found that the 
CME velocity becomes higher after correcting the projection effect. Combined radio, in situ, and white- 
light observations of CME/shocks have also revealed tha t the plane-of-the-sk y speeds for fast halo CMEs 
are always less than or equal to the initial radial speeds (iReiner et all 120071) . 

5 SUMMARY 

In this paper we first analyzed the relation between the white-light intensity and the propagation velocity 
of 3 1 limb CMEs observed clearly by the SOHO/LASCO coronagraph. It is confirmed that slower CMEs 
tend to be weaker in the white-light brightness, meaning that they carry less plasma. We also studied 
the normalized intensity evolution of the CME leading edge along with the radial distance r and fitted 
it with a power-law function, /„ ~ r"'' . It is found that the power index k is averaged around 3.85. 
Considering that the incident emission from the photosphere decreases as r~^, it means that the plasma 
density of a CME frontal loop, n, decreases with r roughly as r^^-^^. Although the cross-section of 
Thomson scattering increases as a limb CME is observed as a halo CME, the abrupt decreases of the 
incident light from the photosphere and n with r lead to that the white-light intensity of halo CMEs 
would be very weak since they should travel a longer distance to be observed by coronagraphs. 

As a further quantitative investigation to resolve the high velocity puzzle for halo CMEs, we per- 
formed a Monte Carlo simulation of 20 000 artificial CMEs with velocity and angular width distributions 
derived from limb CMEs. The Thomson-scattering intensity is calculated for both limb and halo events. 
The simulation indicates that if the limb CMEs with angular width between 20° and 120° are observed 
as full halo events propagating along the Sun-Earth line, a majority of the slower events become so 
weak that their intensity is comparable to the fluctuation of the background solar wind, implying that 
they would fail to be identified by coronagraphs. The average velocit y of the "detectable" h alo CMEs is 
~922 km s^^, which agrees perfectly with previous statistical result ( Michalek et al.ll2003b . We believe 



that the missing of many slow halo CMEs can well explain the high velocity puzzle of halo CMEs. 
As a proof of statement that some CMEs may be neglected even by the state-of-art coronagraphs, here 
we mention the 2007 De cember 7 even t, which was observed by STEREO (Solar TErrestrial RElations 
Observatory) A satellit e dKaiseil 1200 5*). but completely missed by STEREO B since the event was more 
face-on to STEREO B (iMa et al.L 20091) . Statistical work using multi-directional data is of help to clarify 



how many faint halo CMEs are missed and testify our conclusion. 

It should also be emphasized that the halo CMEs in the simulation are full halo ones that originate at 
the disk center and propagate along the Sun-earth line. For the full halo CMEs that propagate aside from 
the Sun-earth line and partial halo CMEs, the white-light intensity decrease would not be so significant 
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as in Fig.|6l Thus, less portion of such events would be missed. In order to quantitatively compare the 
Monte Carlo simulation with observations, these cases would also be considered in a future work. 
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